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Chiral extrapolation of the X (3872) binding energy is investigated using the mod¬ 
ified Weinberg formulation of chiral effective field theory for the DD* scattering. 
Given its explicit renormalisability, this approach is particularly useful to explore 
the interplay of the long- and short-range DD* forces in the A(3872) from study¬ 
ing the light-quark (pion) mass dependence of its binding energy. In particular, the 
parameter-free leading-order calculation shows that the A-pole disappears for un- 
physically large pion masses. On the other hand, without contradicting the naive 
dimensional analysis, the higher-order pion-mass-dependent contact interaction can 
change the slope of the binding energy at the physical point yielding the opposite 
scenario of a stronger bound X at pion masses larger than its physical value. An 
important role of the pion dynamics and of the 3-body DDtt effects for chiral extrap¬ 
olations of the A-pole is emphasised. The results of the present study should be of 
practical value for the lattice simulations since they provide a non-trivial connection 
between lattice points at unphysical pion masses and the physical world. 


PACS numbers: 14.40.Rt, 11.55.Bq, 13.75.Lb, 12.38.Lg 


2 


I. INTRODUCTION 


After more than a decade from the discovery by the Belle Collaboration of the X(3872) 
charmoninm-like state l| its natnre still remains an open qnestion — see Ref. for a review. 
According to the Particle Data Gronp, this state has the mass Mx=(3871.68 ± 0.17) MeV 
and thus resides very close to the neutral DD* threshold with 

Eb = Mdo + Mb*o - Mx = (0.12 ± 0.26) MeV. (1) 

It is therefore natural to assume that its wave function has a lar ge m olecular admixture — 


15l | to mention some, also 


see a vast literature on hadronic molecules, for example. Refs, 
in the context of the X(3872). 

The quantum numbers of the X determined recently by the LHCb Collaboration Q, 
17| | are consistent with its interpretation as an R-wave D^D*^ bound state^ — see, for 
example. Refs. 13l-ll5|. The small binding energy relative to the D^D*^ threshold allows for 
an effective held theory (EFT) formulation of the problem in analogy to the deuteron. The 
pionless EFT based on pure contact DD* interactions was hrst applied to the X(3872) in 


e implications of the heavy quark and heavy havour symmetries were utilised 


20| to predict partner states of the X(3872). However in presence of other 


Ref. Q whi 
in Refs. [3, 

relevant dynamical scales such a treatment is expected to be valid only in a very narrow 
energy region around the threshold. In particular, the 3-body neutral channel D^D^tt^ opens 
at approximately 7 MeV below the D^D*^ threshold while the charged 3-body thresholds 
D^D^tt^ and reside about 2 MeV above it. In addition, the charged two-body 

threshold D^D*"^ is located around 8 MeV above the neutral one. The mass difference 
between the charged and neutral DD* thresholds was shown in Ref. to play a crucial 
role for understanding of isospin violation in the decays of the X into 7r+7r“J/'0 2^ and 

TT~J/'ip 231], for which approximately equal branching fractions were observed. To 
incorporate the long-range pion physics the so-called X-EFT approach was developed in 


Ref. 


24| based on the assumption that pions can be treated perturbatively. Recently this 


ligher-order corrections and then used to predict the 


framework was extended to include 
pion-mass dependence of the X-pole 
energy 26|. On the other hand, perturbative treatment of pions has a smaller range of 


25l | and the hnite volume corrections to the X binding 


^ A proper C-parity eigenstate is always meant by this (and similar) shorthand notation. 
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validity compared to nonperturbative approaches and it has to be used with caution—for 
example, the perturbative framework is known to be not applicable in the deuteron channel 
271] which demonstrates certain similarities with the X(3872). 

The frameworks with nonperturbative pious were employed in many phenomenological 
studies — see, for example. Refs. 28l-l3n| to mention some — all of them, however, include 
one-pion-exchange (OPE) in the static limit, that is under the assumption that the D-mesons 
are inhnitely heavy. Meanwhile, a close proximity of the and 

thresholds to the X(3872) pole suggests that 3-body scales can play an important role in 
this state, so that neglecting the 3-body dynamics one distorts the analytical structure of 
the amplitude in the kinematical region of interest. It has to be noticed, however, that the 
proper inclusion of the 3-body dynamics requires special care. For example, it is shown in 


Ref. 


3l[ | that the 3-body unitary cuts play very important role in the D^Dis system, if the 


Djs width is dominated by the R-wave —)■ decay. In particular, it is demonstrated 

that if the D/? —)■ Dan coupling is sufficiently strong to produce a bound state |32|, ISSj, it 
is, at the same time, necessarily sufficiently strong to provide the state with such a large 
width that it becomes unobservable. In turn, in the case of P-wave vertices, the system 
at hand demonstrates additional difficulties since the selfenergy loops diverge, so that the 
system requires a proper treatment to avoid false conclusions — see an example of such 


conclusions in Ref . 1341 and its detailed discussion in Ref. [3^. In particular, contrary to 
the claims of Ref. 3J|, it is shown in Ref. 35| that the one-pion exchange (OPE) potential 
in the DD* system is well dehned in the sense of an effective held theory only in connection 
with a contant operator. 

In Refs. 


36l. l37l |. the properties of the X(3872) molecular state were studied in a heavy- 
meson EFT framework with nonperturbative pions including all relevant 3-body scales. It 
was understood that the dynamical treatment of pions had a big impact on the X line 
shape and, in particular, on the partial decay width X —)■ Furthermore, it was 

shown in Ref. 36(| that the static OPE approximation was not adequate to analyse the role 
of the long-range pion dynamics in the X(3872), since it corresponded to an uncontrolled 
modihcation of the proper dynamical scales related to the DDn cuts and to the neglect of 
the imaginary part of the DD* potential. Meanwhile, the role of nonperturbative effects for 
these observables appeared to be quite moderate, as follows from the agreement between 

n nn 

the results of the nonperturbative calculations [36| and those in the X-EFT [l4 l24l|. 
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In Ref. 


37l |. the nonperturbative framework developed in Ref. 36|] was generalised to 


study the dependence of the X binding energy on the light-quark mass or, equivalently, 
on the pion mass. The use of non-perturbative one-pion exchange for chiral extrapolations 
allows one to extend the region of applicability of the approach to larger pion masses which 
is important for analysing the results of lattice QCD calculations. 

In this work we address another important issue which is related to nonperturbative 
renormalisation of the 3-body Lippmann-Schwinger or Faddeev-type equations to describe 
the interaction between heavy mesons in the X. The standard nonrelativistic approach 
to heavy mesons leads to coupled-channel integral equations for the scattering amplitudes 
which, at leading order in the EFT expansion, are linearly divergent. As a consequence, it¬ 
erations of the truncated potential within the dynamical equation generate an inhnite series 
of ultraviolet (UV) divergent higher-order contributions to the amplitude which cannot be 
absorbed into a hnite number of counter terms (contact interactions) included in the poten¬ 
tial. In other words, the coefficients in front of the logarithmic and power-law divergences 
appearing in the iterations of the equation involve powers of external momenta which can 
only be removed if an inhnite number of higher-order (derivative) contact interactions is 
included. The standard way to deal with this problem is to employ a hnite UV cut-oh of 
the order of a natural hard scale in the problem which would suppress the unwanted higher- 


order contributions, as advocated in Ref. 


This strategy was followed, in particular, in 


Refs. [ 19 , 20, 36, 37, 391. Exactly the same problem with renormalisation emerges also in 
the context of nuclear chiral EFT — see, for example. Refs. 40|, l4]| and references therein. 
In particular, a hnite cut-oh was employed for the construction of the NN potential and 
the few-body nuclear forces within chiral EFT — see Ref. 42| for a review. This procedure 


induces cut-oh artefacts which might turn to a nontrivial issue, in particular, for chiral ex¬ 
trapolations since it might be difficult to control the pion mass dependence of short-range 
interactions in a systematic way. Note also that in Ref. the m 7 r-dependence of the con¬ 
tact interaction was promoted to the leading order to maintain the renormalisability of the 
scattering amplitude at unphysical pion masses. 

Recently a novel, renormalisable (in the EFT sense) approach to nucleon-nucleon scat¬ 


tering with nonperturbative pions was proposed in Ref. [43|. S 


invariant form of the ehective Lagrangian, the authors of Ref. 


arting from the Lorentz 


43| derived a 3-dimensional 


dynamical equation which complies with the relativistic elastic unitarity and which is renor- 
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malisable at leading order of EFT. Indeed, in the snggested approach, all logarithmically 
divergent contributions generated by iterations of the potential can be fully absorbed into the 
redehnition of the leading-order contact terms. Then higher-order contributions are subject 
to a perturbative treatment in this approach. It should be stressed that the central point of 
the approach is noncommutativity of the nonrelativistic expansion and the renormalisation 
procedure while after renormalisation relativistic effects as such provide only minor impact 
on the low-energy observables, as it should be in EFT. Apart from its transparency with 
regard to renormalisation, this approach allows one to remove hnite cut-off artefacts and it 
is very well suited for carrying out chiral extrapolations and studying correlations between 
the effective ran ge p arameters induced by the analytic structure of the long-range forces 


— see Refs. 4J, l45| for the corresponding results in the NN sector. Given the same UV 
behaviour of the dynamical equations for NN and heavy meson-antimeson scattering, we 
demonstrate that the method of Ref. ^ can be used to reformulate the nonrelativistic 3- 


body approach of Refs. 3a, |37|] in terms of renormalisable integral equations. We apply the 


resulting theoretical framework to study the quark mass dependence of the X(3872) bind¬ 
ing energy. In contrast to the hnite cut-off formulation, the mjr-dependence of the binding 
energy is predicted at leading order in a renormalisable approach. 

The paper is organised as follows. In Sec. [TTl we give a brief introduction to the method 
suggested in Ref. 43|]. In Sec. IIIII the problem of the DD* interaction is formulated in a 
closed selfconsistent form which makes it possible to appeal to the approach discussed in 
Sec. m In Sec. m we present and discuss the results of our calculations. We summarise our 
hndings in Sec. |Vl All necessary technical details are collected in Appendix |A1 


II. NON-RELATIVISTIC (LIPPMANN-SCHWINGER) EQUATIONS VERSUS 
EQUATIONS WITH RELATIVISTIC UNITARITY 


A proper nonrelativistic expansion of low-energy physical quantities can be done by cal¬ 
culating these quantities in a Lorenz-invariant theory and expanding the hnal result in the 
powers of the velocity y (we work in the natural system of units setting the speed of light 
c = 1), see e.g. Ref. [4^ where this issue is discussed in the one-nucleon sector. On the other 
hand, one can perform the nonrelativistic expansion at the level of the Lagrangian of the 
theory. However, this expansion does not commute with the loop integration. This can be 
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exemplified by a simple calculation adapted from Ref. . Consider a scalar two-point loop 
function which is logarithmically divergent and therefore should be regularised. With the 
simplest regularisation prescription given by a sharp cut-off in the 3-dimensional momentum 
it reads 




/ d^ke{k-\k\) 

(27r)'^ J [k^ — m? + iO] [{P — ky — m? -\- zO] ’ 

where P = -|- p^, 0 ). The integral can be evaluated analytically for A > |p| with the 

result 


I = 


t\p\ 


+ 


\P\ 




: In 


+ p"^ + |p|\/A2 -|- m? 


m 


k? — p"^ 


1 A -|- y k^ + m? 

— In-. 

TT"' m 


( 3 ) 


A nonrelativistic expansion of the integrand in Eq. ([2]) implies the strong inequality 
IpI A <C m that is equivalent to the 1/m expansion made prior to the 1/A expansion in 
the exact result ([3]), that yields 

A 


/ = 


i\p\ 

27rm 


TT^m 


-I- ..., 


( 4 ) 


where the ellipsis denotes suppressed terms. Divergence in expression (jl]) is linear, that is, it 
is stronger than that in the original integral (|2]) which is a consequence of the nonrelativistic 
expansion of the integrand. 

On the contrary, keeping the integrand relativistic and performing the nonrelativistic 
expansion after integration is equivalent to imposing a different, and more natural, strong 
inequality |p| -C m -C A and, therefore, the 1/A expansion is to be performed in Eq. ([3]) 
before the 1/m expansion. This leads to a different result for the real part of the integral, 

i\p 


r 1 , 2A 

/ — —-^ In-h . . . , 

27rm TT^ m 


( 5 ) 


which reveals the logarithmic divergence, in agreement with the UV behaviour of the original 
integral. 

Thus nonrelativistic expansion of the integrand changes its ultraviolet behaviour and the 
final result differs from the relativistic expansion of the exact expression for the integral. 
This difference is caused by noncommutativity of the nonrelativistic expansion and the 
loop integration. Because of this noncommutativity, in order to reproduce the results of 
the Lorentz-invariant theory, one needs to add compensating terms to the nonrelativistic 
effective Lagrangian. Therefore, more singular behaviour of the nonrelativistic equation 
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leads to perturbative nonrenormalisability already for the leading order (LO) potential. In 
particular, iterations in the Lippmann-Schwinger equation generate power-law divergences 
with coefficients of a progressively increasing power of the momentum/energy. As the LO 
potential does not contain momentum-dependent contact interactions, one cannot get rid of 
these divergences by absorbing them into redehnition of the parameters of the LO potential. 
Adding any hnite number of momentum/energy dependent terms does not resolve the issue. 
While this is not a problem when calculating a hnite number of diagrams, it is rather 
disturbing when solving integral equations. Except for some trivial cases, it is not possible 
to take into account contributions of an inhnite number of compensating terms required 
for “correcting” an inhnite number of iterations. One is, therefore, forced either to keep 
the ultraviolet cut-oh hnite (A ~ m) or to resort back to the original Lorentz-invariant 
formulation of the theory, although the ehect of relativistic corrections at low energies is, of 
course, small after renormalisation. 

On the other hand, iterations of the Lippmann-Schwinger equation without nonrelativistic 
expansion generate only logarithmic divergences. This guarantees a perturbative renormal- 
isability of the theory at LO, that is, all divergences can be removed by renormalising the 
coupling constant of the LO contact interaction. 

In Ref. 37|, the nonrelativistic Lippmann-Schwinger equation was solved that corre¬ 


sponds to the nonrelativistic expansion of the integrand as was explained above. In the 
present paper we deal with the relativised Lippmann-Schwinger equation and, therefore, the 
LO amplitude is obtained by solving a renormalisable integral equation. Analogously to the 
nucleon-nucleon scattering in the modihed Weinberg approach of Ref. 43|, the relativised 
integral equation for the DD* system has a milder ultraviolet behaviour if compared to 
the nonrelativistic Lippmann-Schwinger equation. It has to be noticed, however, that the 
integral equation becomes nonrenormalisable if higher-order corrections to the leading-order 
potential are also treated nonperturbatively. In particular, by iterating higher-order con¬ 
tributions in the potential one generates divergences with such structures of momentum- 
and/or energy-dependent coefficients which are not present in the iterated potential, that is 
these divergent contributions cannot be absorbed into the redehnition of the contact terms 
included in the potential at the given order. On the other hand, renormalisability is retained 
by treating corrections perturbatively. In particular, if we denote the LO amplitude as Tq 
and the NLO corrections to it as Ti, that is we have the following perturbative expansion 




of the full amplitude: 


r = r„ + eri + 0(£"), 


( 6 ) 


then the inverse amplitude takes the form 

T-i = T,-\To -eT,)T,-^ + 0{e%-^), (7) 

where £ stands either for the expansion parameter of chiral EFT, e: ~ with 

A^ being the chiral symmetry breaking scale, or it corresponds to the expansion around the 
physical pion mass that is £ ~ (m^ — In what follows, while we stick to 

the leading-order chiral potential, we investigate the pion-mass dependence of the X(3872) 
binding energy including corrections at the NLO which appear as one goes away from the 
physical point. 

Note that the expression of Eq. (iTl) gives an explicitly unitary amplitude, however it also 
includes selectively resumed higher-order contributions which do not affect renormalisability 
of the scattering amplitude. A bound state corresponds to the zero of the inverse amplitude 
([7]). Finally, we use the superscript “ph” to label quantities taken at the physical point, that 
is for ttItt = 


III. SYSTEM OF COUPLED-CHANNEL INTEGRAL EQUATIONS FOR THE 

DD* PROBLEM 


In this Section we outline briefly our theoretical formulation of the problem. We follow 
the lines of Ref. [37| adapting the approach according to Ref. [43| . 

The lowest-order D*Dn interaction Lagrangian is taken in the form 241 


C = 


I■ VttV^R) + V7r“ ■ 




The dimensionless coupling constant Qc is related to the D*^ decay width as 




r{D*° = 




24:71 ’ 


( 8 ) 

(9) 


where q = A^/^(m^Q,mQ,m^o)/(2m*) is the center-of-mass 3-momentum of the outgoing 
particles and X{x,y,z) is the standard triangle function — see the dehnition in Eq. flA6p . 
Here and in what follows, m*, m, and denote the masses of the D* meson, D meson. 
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and pion, respectively. Charged and nentral states are distingnished by an additional index, 
for example versns m*o- 

The DD* potential at LO in Chiral Effective Field Theory (ChEFT) consists of the OPE 
and the S'-wave derivativeless contact interaction Co, 


p') = (p + pT(p + p')”'f’«(p. p') + Co s""'. 


( 10 ) 


where indices n and n' are contracted with the corresponding indices of the D* polarisation 
vectors. Here 

and Dgl] and {i,j = 0,c) stand for the DDtt and D*D*7r propagators written in the 
framework of the Time-Ordered Pertnrbation Theory (TOPT) — see Fig. [H 


E^o{p + p') (Eoi(p) + Eoi{p’) + E^o{p + p') - m\ i = j, 

DSil{P,p')={ ( 12 ) 

En-{p + p') [Eoiip) + Ed^ ip') + E^cijp + p')- My j, 


E^o {p + p') (Ed* (p) + Ed* (pO + E^o {p + p') - m\ i= j, 

DS]iP,P’)={ ^ ^ ( 13 ) 

E^c{p + p') (Ed* (p) + Ed* ip') + E^c{p + p') - My j. 

For convenience, the energy E is connted relative to the nentral two-body threshold. 


M = m*o + mQ + E, 


(14) 


while the energies of the individnal particles are 


EnAP) = yp^ + EdAp) = yp‘^ + ml Ed*{p) = yp^ + ml. (15) 

The OPE potential flTU]) interrelates the fonr H-meson channels dehned as 

|0) = |0) = D^D*^, |c) = D+D*-, |c) = T)~T)*+. (16) 

Then the system of conpled-channel Lippmann-Schwinger eqnations for the DD* t-matrix 
elements a^AiPiP') (^^oAPiP') C-even channel has the form 


aSo”'(P.P') = W(P,P') I ar'(fe,P'), 


i=0,c 


^cO 


'(P.P') = KVp{p,p') - A. / d^kVr{p,k)y^aZ'''(k,p'), 

i=o,c 2 d ) 


(17) 
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FIG. 1: Diagrams in time-ordered perturbation theory corresponding to the (inversed) 3-body 
propagators (left plot) and (right plot) {i,j = 0,c). The double and single solid lines 
refer to the D* and D, respectively, while the dashed lines refer to pions. The thin vertical line 
pinpoints the intermediate state. 


where Aq = (Olri • r 2 |()) = (clri • r 2 |c) = 1 and Ac = (0|ri ■ r 2 |c) = (cIti • r 2 |()) = 2 are the 
isospin factors for the vr^- and 7r^-exchange, respectively. 

The partial wave projections of potential (1T0|1 on the relevant ^Si and ^Di partial waves 
read {x = cos 6 where 6 is the angle between the momenta p and p') 


P') = Co + - I Fij{p, p', x) + 2pp'x) dx, 




rSDr _ l\ 


vrip.p') = - 


6 7-1 

V2 r 


Fij{P:P\x) 

Fij{p,p',x) 




v^rip.p') = - g 


Vif^{p,p') = ^ j^Fij{p,p',x) 


p'^ + ~ ^ ^ 

P^ + p'"^ ( ~ ^ ^ 2pp'x 


dx, 

dx, 


/ 2 f2^ 1 ^ 2 11 / 9 f 3 3 

[p +p )(-x - 2 ) + Yo^P ^ + Yo^P ( 2 ^ ~ 2 ^ 


dx. 


Because of the P-wave nature of the D* —)■ Du vertex, the Dir loop operator 
S(s, m*, m, tUtt) diverges and it is subject to renormalisation. The necessary details of the 
renormalisation procedure are given in Appendix El while here we only quote the final result 
for the inverse two-body propagators Aq and Ac entering system of equations flTTl) : 

Ao(p) = 


Ac(p) = 


P^o (p)Po*o(p) 

Edo{p) + Ej:,*o{p) - M 


[ c 

Sij(s, m*o, m^o, mo) + 2Er(s, m*o, 


ijp) 

Edc{p)Ed*c{p) 

Ed<={p) + Ed*<={p) - M 

rricm^c 

[ C 

Sij(s, m*c, m^o, me) -b 2T.r{s, m*c, m^^, mo) 


1*0 


2Eo..[p) 


( 18 ) 






























11 


where 


3 + 2Ed*{p){M - Eo{p) - En-ip)) 


( 19 ) 


for the off-shell D* resonance and ^r{s, m*, m, m^) is the renormalised loop operator defined 
at the “running” pion mass, that, in particular, brings about the quantity (, 


r' = 1 - 


9r 


In 


mz 


3847r2 (^f)2’ 

with the renormalised coupling constant dehned as (see Appendix 

y/mom^o 


( 20 ) 


Qr — Qc- 


f. 


( 21 ) 


For future discussion of the static approximation, we also consider a simplihed case of 
the constant width which corresponds to the substitution of the constant s = ml in Eq. flTSl) 
instead of the “running” s, as given in Eq. (IT^ . 

We are now in a position to introduce the power counting in the parameter ^ 

The tUtt- dependence of the coupling constant gc is extracted from Ref. [48j and is discussed 


in Ref. 37|]. In particular, at LO Qc remains constant while at NLO it acquires corrections 
of the order of m^. Similarly, /^, the masses of the D and D* mesons, and the decay width 
r{D*^ —>■ D^'^) take their respective physical values in the LO calculation. The central 
issue of this work is related to the m,r-dependence of the contact interaction Cq — see 
Eq. (ITOll . Since the nature of this interaction is obscure, the dependence can only 

be guessed using the principle of naturalness. Below, we discuss in detail generalisation of 
the corresponding approach developed in Ref. Meanwhile, regardless of the particular 
source of the dependence C'o(m.n-) it only appears at NLO, so that the contact interaction 
remains constant at LO and, therefore, the problem is fully hxed to provide a prediction for 
the behaviour of the X-pole as the pion mass leaves the physical point. Furthermore, the 
pion-mass dependence at LO occurs only due to the pion energies in the DDti propagator 
and pion-mass effects in the renormalised loop — both are the equally important parts of 
the three-body DDti dynamics. To hnalise the setup of the problem, we quote the values of 
various parameters at the physical pion mass used in the calculations. In particular, in the 
physical limit of one has = 92.4 MeV, then 

m^o = 134.98 MeV, mf = 1864.84 MeV, = 2006.97 MeV, 

= 139.57 MeV, mf = 1869.62 MeV, mf = 2010.27 MeV, 
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FIG. 2: Pion mass dependence of the 2f(3872) binding energy at LO. The results of the full 
dynamical theory with 3-body effects included (black solid curve) are confronted with the simplified 
formulation with static OPE (blue dotted curve). 

and the values 


^ph(^*o ^ ^ 42 keV, rP^(D*° ^ L)° 7 ) = 

can be deduced from the data for the charged D* decay modes 


the couplings introduced above are 


37| 


21 keV (22) 

The physical values of 


gf = 0.61, gf = 12.7. 


(23) 


IV. RESULTS AND DISCUSSIONS 

We are now in a position to discuss the results for the pion mass (or, equivalently, light- 
quark mass) dependence of the X(3872) binding energy Esim.,,). We start from the discus¬ 
sion of the LO results. As was explained above, the contact term Cq is m^-independent at 
this order, so once it is adjusted to reproduce the binding energy at the physical pion mass 
(for dehniteness we set = 0.5 MeV), the scattering amplitude can be calculated 

for unphysical pion masses without loss of renormalisability of the LO equations. Therefore, 
at LO of our EFT, the dependence Es(m 7 r) can be predicted in a parameter-free way — 
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see Fig. [2l At this order, the pion mass dependence of Eb originates only from the pionic 
effects in the OPE potential and from those in the renormalised selfenergy loops — see 
Eq. flTSll . The binding energy at LO demonstrates a clear tendency to decrease with the 
growth. Note that a similar behaviour of the binding energy was observed in Ref. 441 
for the deuteron. Furthermore, the slope of the binding energy in at the physical point, 
exhibits a strong sensitivity to the 3-body DDtt effects. In particular, 
neglecting the 3-body dynamics (the so-called static OPE) results in a much steeper fall of 
the binding energy — compare the dotted (blue) line versus the solid (black) in Fig. [2l 
Since no real experiment is possible for unphysical pion masses, the only source of in¬ 
formation on the X pole fate for the m^’s exceeding the physical pion mass is provided by 
lattice simulations. Such calculations are indeed being performed and most of them predict 
an increase of the binding energy with the growth. For example, different lattice collabo¬ 


rations observe this type of behaviour for the deuteron — see, for example. Refs. 49|-^ and 
references therein^. Also the hrst lattice calculations for the X(3872) indicate the existence 
of a stronger bound X for 53l-l55l|. Although these results still suffer from poten¬ 

tially large hnite-range corrections, as pointed out in Ref. 2^ within X-EFT, they raise an 
important question of whether such a behaviour of the binding energy can be understood 
theoretically. To this end, we go beyond LO and proceed to NLO thus including corrections 
quadratic in In particular, we allow for an m 7 r-dependence of the short-range interac¬ 
tion which therefore goes away from its physical value. Thus, we consider (for simplicity all 
indices are omitted) 


Rnlo = Vope{Pt p\ 0 + Co + D{^‘^ - 1), ^ 


(24) 


where the hrst two terms on the right-hand-side stand for the LO potential (ITOll while the 
last term accounts for our ignorance of other dynamical scales but those related to the 
OPE. As was discussed in the previous Section, renormalisability of the theory requires 


all operators beyond the LO to be included perturbatively. Following Ref. 


37|. we hx the 


unknown coefficient D to the slope of the binding energy at the physical pion mass which 
is therefore considered as an additional input quantity. For example, in Fig. [3] we illustrate 
the behaviour of the binding energy for the slope {SEb/ dmD\^^^^ph ^ E "^— see 


On the other hand, the HAL QCD Collaboration found no bound state in the NN ^Si-^Di channel 52l |. 


2 
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the dashed curve in the left panel. While the sign of the slope was hxed to provide a 
growth of the binding energy with the pion mass, its magnitude was chosen to comply with 
naturalness. Specihcally, we assume that the shift of the binding energy 6Eb Ef for 
SuIt, ~ can be interpreted as natural. Indeed, the slope predicted at LO due to OPE 
fulhls this criterion: {OEb/ dmD\^^^^ph ^ —1.5E'^/m^^‘^. Therefore, to study the case of 
a stronger bound X, we £x the slope to be {dEB/dml)\^ ^ E'^ 1171^“^. Interestingly, 
in a theory with the same polynomial behaviour of the contact operator but without pions, 
one would observe a much flatter behaviour for the same slope {dEB/dm^)\^ 

as shown by the dashed-dotted curve. The difference between the two curves demonstrates 
the role of dynamical pions as an explicit long-range degree of freedom. As seen from 
Fig.[3l the contact interaction provides a smooth background for a rapidly varying pion-mass 
dependence stemming from OPE. Therefore, integrating out pions and the corresponding 
3-body soft scales while still trying to, at least partially, compensate for neglecting these 
long-range effects, one would inevitably arrive at unnaturally large m.n--dependent coefficients 
accompanying short-range operators. 

On the other hand, one may question a justihcation of the perturbative inclusion of the 
m^-dependent short-range interaction in Eq. fl241) . Given the shallowness of the physical 
X state, even a small variation of the slope within its natural range, as discussed above, 
has a sizeable impact on the dependence of the binding energy. In order to verify the 
validity of the perturbative approach, we employ resonance saturation to model higher-order 
contact interactions by means of a heavy-meson exchange. In particular, we consider the 
NLO potential in the form 


Gnlo — hbpE(p, 0 + C'o + [3- 


9c 


[£-q){e*-q) 


(47r/^)2 q^ + M^ + 5M^ - 1) ’ 

where q = p + p' and the scale M is varied in the range M = 600...800 MeV that cor¬ 
responds to a typical heavy-meson mass. The parameter /3 accounts for the difference in 
the strength of the heavy-meson exchange potential relative to that of OPE. It is expected 
to take values around unity and it could be, in principle, adjusted to the DD* effective 
range. However, given that the latter is unknown, f3 is varied within a suitable range of 
values from 1 to 2 which we treat as natural. The term — 1) in the denominator 

accounts for the pion mass dependence of the heavy meson mass^ with 6M adjusted to 


Notice, a similar EFT approach based on the resonance saturation hypothesis was used in Ref. 56| to 
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the slope {OEb/ dm\)\^ =TnF^- have verihed that for 5M = 0, the dependence EsirriT^) 
is basically indistinguishable from the LO one that conhrms the results to be insensitive 
to the details of the short-range interaction, as expected. The form of the NLO potential 
(j25|l ensures that the corresponding scattering amplitude is renormalisable so that the NLO 
calculations can be carried out in the same way as the LO ones. Then, hxing the slope as 
before, {dEB/dml)\^^^^,j. ^ E^/m^'^, we obtain the dotted (red) band in Fig. [3l A very 
good agreement between this band and the dashed curve in the considered range of pion 
masses conhrms that higher-order effects originating from the non-perturbative resumma¬ 
tion of pion-mass-dependent short-range interactions are minor. Therefore, the perturbative 
treatment of the pion-mass-dependent short-range interaction, as given by Eq. fl2Tll . is in¬ 
deed justihed even though it brings about a new effect — the slope of the binding energy 
may change its sign compared to the LO result depicted in Fig. [2l 

It is also instructive to compare the results of the present study with those obtained in 
Ref. 371] in the heavy-meson formulation with a hnite cut-off — see the right panel in Fig. [31 
At NLO, both approaches are consistent with each other as may be expected since both 
formulations are justihed in general. Meanwhile, as was already discussed above, the LO 
equation is explicitly renormalisable and predictive in the present formulation in contrast to 
the approach of Ref. ^ , where the m^-dependent contact interaction requires an additional 
input to be included to maintain the renormalisability of the scattering amplitude. In 
addition, the requirements of naturalness are much easier to formulate and to apply in the 
current approach since the dependence of the results on the cut-off is eliminated. Thus, we 
emphasise that the role of pion dynamics can be understood in a much more transparent 
way using the explicitly renormalisable theory which is free of hnite cut-oh artefacts. 

To further clarify the role of the dynamical pions at NLO, we assume that there exist 
gedanken lattice data at unphysical pion masses. These data could be used to adjust the 
parameters of the short-range potential. Then, once the short-range physics is hxed, the 
theory can be extrapolated to the physical point in and confronted with the experimental 
data. In particular, if the lattice calculations provide two measurements of the binding 
energy of the X made for two unphysically large pion masses then the suggested approach 
allows us to establish the correct extrapolating formula to the physical point and thus to 


constrain the pion mass dependence of the short-range NN forces. 





1 R 



FIG. 3: Pion mass dependence of the X(3872) binding energy at NLO. Left panel: the dashed 
line is for the perturbative treatment of the rriTr-dependent contact operator at NLO while the 
red dotted band represents the nonperturbative results employing resonance saturation. Right 
panel: comparison of the results obtained in the heavy-meson formulation of Ref. 371] with the 


finite cut-off A E [500 MeV, 700 MeV] (black hatched band) with the nonperturbative results of 
the current study employing resonance saturation and the cut-off A —>• oo (red dotted band). The 
dashed-dotted line in both panels corresponds to the calculation without pions. The (blue) dot 


with the error bars shows the result of the lattice calculation of Ref. 




predict the corresponding valne EB{rn'^). In addition, information on the behaviour of the 
short-range interactions in the X, which can, in this way, be extracted from the lattice 
data, may shed light on the nature of the binding mechanisms in the X. This establishes 
an important link between the EFT approach and lattice simulations for hadronic molecule 
states. 

As an illustration, the chiral extrapolations for the two theories, the one with dynamical 
pions and the one with the static OPE, are compared with each other in Fig. 01 For 
definiteness the gedanken lattice result is taken at m-„ = as indicated by the arrow 

in Fig. m If the slope is chosen such that the theory with dynamical pions provides the 
correct extrapolation to the physical pion mass, the static theory with the same slope yields 
a signihcant overbinding of the X(3872), by more than a factor of 3 in the binding energy. 
In addition, one can see from Fig. 0]that the extrapolation curve from an unphysically large 
pion mass, close to the values used on the lattices, to the physical point is nontrivial and the 
corresponding extrapolating behaviour has a strong curvature. This illustration emphasises 
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FIG. 4: Pion mass dependence of the 77(3872) binding energy. The red dotted band is for the full 
calculation with dynamical pions at NLO while the blue crossed band is for the static OPE. 

the importance of the 3-body effects for the chiral extrapolations for the X and pnts in 
question possibility of using any simple ansatz for the extrapolation formula. 


V. SUMMARY AND CONCLUSIONS 


In this work we developed an explicitly renormalisable framework to study chiral extrap¬ 
olations of the binding energy of the X(3872) beyond the physical pion mass. This approach 
is free of the hnite cut-off artefacts which is a precondition for a systematic control over the 
pion-mass dependence from the short-range interactions. The pertinent results of our work 
can be summarised as follows. First, the interplay between the long- and short-range forces 
in the X appears to be quite nontrivial, as was already pointed out in Ref. 37|]. If the X turns 
out to be less bound for the pion masses exceeding its physical value, the m^r-dependence 
of the X binding energy is entirely governed by the explicit pion-mass dependence of the 
OPE potential. On the other hand, a stronger bound X would signal the importance of the 
mTT-dependent short-range interactions in addition to pionic effects. Confronting our results 
with those of the lattice simulations could allow one to extract valuable information on such 
short-range interactions and, possibly, to disclose the nature of the binding forces in the X. 

Secondly, our Endings are of a practical value for the lattice simulations since they open 
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the way to override the gap between the unphysically large pion masses used on the lattices 
and the physical limit. It follows from our results that the corresponding interpolating curve 
has a strong curvature and it is strongly affected by the 3-body effects in the X. 

Last but not least, the approach developed in this paper can also be adapted to other 
near-threshold states, the X(3872) being just the most prominent and therefore the most 
extensively studied one. 
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Appendix A: The Dtt loop operator and the D* propagator 


Gonsider an unstable vector mesonic state which decays in the P wave into a pair of 
(pseudo)scalar mesons. We start from its inverse propagator 

T)(s) = s - mgS(s), (Al) 

where s is the invariant energy (s = p^), rrio stands for the ”bare“ mass, and E(s) denotes 
the selfenergy loop operator. The one-loop contribution to the selfenergy has the form'^ 

ns) = -gll,ue>^e*'', (A2) 


where Qq is a dimensionless bare coupling constant and e’s stand for the polarisation vectors 
of the unstable meson such that {e ■ e*) = —1 and {e ■ p) = 0. 

The loop integral corresponding to the decay of the vector meson into two mesons of the 
masses mi and m 2 reads 




d^k 


knkj. 


(27r)’^ {IP — mf)((fc -|- p)2 — m^ 


— + B 


Pl^Pv 


pZ 


(A3) 


4 


Vector meson self-energies in chiral EFT were first discussed in Ref. 
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and, due to the property (e-p) = 0, only the function A is relevant for the loop operator 
Since the loop integral flA3|) diverges quadratically in the limit of n —)■ 4, the coefficients Aq 
and Ai in the Taylor series of the function A(s), 


A(s) — Aq + AiS + Areg(s), 


(A4) 


are singular in this limit, while the residual function Areg(s) is regular. By a straightforward 
calculation in the dimensional regularisation scheme, it is easy to hnd that 


An — 


1927r2 


((2m2 + m^)C(m2) + (2771^ + m2)C'(mi)) , Ai = 


1 


3847r2 




12s 


where 


/2®®(s,mi,m2) = 


771^ 77I2 , TTli / 2 2'i rreg/ \ 

2 In-h A(s, m^, 1712)12 s) 

lOTT^ 7712 


m? — TTlo , 777.1 

-1 H- - -^In — 


(C(777i) + C(7772)) , 

(A5) 




7772 


.^72,2) ^ _|_ ^2 _ g _ yl/2(^g^ 777|) 

S 2777 i 7772 


and 


/1 \ 77T^ 1 

C(777) = f - - log(47r) + 7 i 5 - Ij +ln—, e =-(4 - 77 )0, 7 ^; = -r'(l) cs 0.5772, 

while /i is the scale in dimensional regularisation and the triangle function is dehned in the 
standard way, as 

A(x, 7/, z) = + 7/^ + — 2x7/ — 2x2; — 2yz. (A6) 

Then, the one-loop contribution to the selfenergy in Eq. (IA2p takes the form 


^(■5) — fi'oA('S) — PoiAn + Ais - 1 - Ai.eg(s)), 


(A7) 


and it is subject to renormalisation which we perform by expanding S(s) near the renor- 
malised vector meson mass 


S(s) = ReS(777^) -|- ReS'(777^)(s — 777^) -|- Sreg(s), 


(A8) 


® In general it is preferable to relate the renormalised mass to the complex pole of the propagator 581163 1. 
However, at the one-loop level it is sufficient to use the real part of the inverse propagator since the 
difference between the two approaches occurs only starting from the two-loop order. 
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where 


ReS(m^) = gl{AQ + Aim\ + ReA,eg{m%)), = gl{Ai + ReA'^^^imj^)), (A9) 

while 

Si.eg(s) = S(s) - ReS(mK) - ReS'(mK)(s - m%) = glAR{s), (AlO) 

with 

Ar{s) = Areg(s) - ReAreg(m^) - ReA;^g(m^)(s - m\). (All) 

Notice that A/j(s) is hnite and it does not depend on the auxiliary regularisation scale /i. 
Dehning the renormalised mass such that — ReS(m|j,), we have: 

D{s) = {s-ml + ReS(m^)) + ReS'(m^)(s - m^) + Sr(s) = Z~\s - m^) + g^ARls), 


with 


Z ‘ = 1 + ReE'(m|) = 1 + g§(Aj + ReA' (m%)), 


(A12) 


where it was used that Ai is real. 

Consider now the combination entering the system of equations for the scattering ampli¬ 
tudes fflTj) . 


9i 


al 


9r 


D{s) Z-^{s-m%) + glAR{s) s - m|^^(s) 
where the renormalised coupling constant gR is dehned as 


> Ejf(s) — gjiAji(s)y (A13) 


al = Zal (A14) 

In particular, for the case of the an obvious identihcation of the parameters is 
niR = m*o, '^1 = 1^0 , 7712 = 771^0. Then, using Eq. ([H]) and the standard relation between 
the loop operator and the width, 

-)■ D°7r°) = -^RuEr^s = mlQ,mR = 171^0,771^0, mo) 

m^o 

= ^ImAreg(s = m^o, mR = m*o, m^o, mo), 
where ImAreg can be found from Eq. flASD to be 

g3 ^ 

ImAreg(s = m%,mR = m*o,m,,o,mo) = 7 --, q = - -mg, m^o), 

Z^TlTTlii^Q ZrUi^Q 


(A15) 


(A16) 
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one arrives at the following relation between the conplings qr and Qc- 

^mom^o 


9r — 9c- 


/, 


(A17) 


that completes renormalisation programme at the physical point. 

Away from the physical value of the pion mass one can write: 

D{s) = s — ml + S(s) = s — m\ + E(s) — ReS(m^) = s — + gl A(s) — ReA(m^) 


= + Ai{s - m|) + Areg(s) 


R,6ylj-0g 


(A18) 


= S - ™R + 9o Mo - M.) + ReA^g(mn)(s - irij,) + i„(s) 


where 


Ak(s) = A„g(s) - ReA„g(m|) - ReA'P^(mJ,)(s - mj), 


(A19) 


that is, Aji{s) is dehned with the derivative in the subtracted term evaluated at the physical 
point. This allows one to avoid m^-dependence in the renormalisation factor Z and to 
preserve its dehnition in the form of Eq. (]A12p . Then 


D{s) = (s - ml) + gl{Ai + ReA[f^{ml))j + glAnis) 

= (s - ml) + gl{Ai - + gliAf" + ReA'Pg(m^)) j + glAR{s) 

= {s — ml)(^gl{Ai — A^'^) + Z + glAji{s) (A20) 


= Z 


-1 


{s-ml){l + gl{Ai- 


+ 9rAr{s) 


= Z 


-1 


2 ^ / 1 9l rn ^ 


(s - m^) 1 - 


+ 9rAr{s) 


where the dehnition of the renormalised coupling gR, Eq. flA14D . and the explicit form of Ai, 
Eq. (]A5p . were used. Therefore for the “running” pion mass instead of Eq. (1A13P one has 


9i 

D{s) 


9l 


9r 


9r 


'-1 


{s - ml)C-^ + glAR{s) 


(s - + E/j(s) Dr{s) 


where 


and 




9l rn 


3847r2 


(m 


ph'i 2 


^h(s) = 9l [74reg(s) - ReAreg(m^) - ReA;P^(m^)(s - ml)] . 


(A21) 


(A22) 


(A23) 
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After integrating over the zeroth component in the loop (one picks the D-meson to be on 
energy shell), one hnds 


■s - rnlQ Ri 2 Ed*o (p) (po - Ed*o (p) ) = 2Ed*o {p){M - E^o (p) - Ej:,*o (p)). ( A24) 


Therefore, one arrives at the formula 


Dd*o (p) 


2Ed,o (p) 

c 


M - Er)o{p) 


Ed*o{p) 


+ ^r{p)j 


(A25) 


for the renormalised inverse D*^ propagator valid for unphysical pion masses. Obviously, a 
similar formula holds for the inverse propagator for the charged particles, as well. 

The inverse two-body propagators Aq and entering system of equations ffT7|) can be 


written as 


Ao(p) 


Epo (p) T>d*o(p) 
mo 2m*o 




Ep-jp) T>d*c(p) 
me 2m*c 


(A26) 
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